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Abstract. We establish an estimate for the fundamental solution of the heat equation on a closed Rie- 
mannian manifold M of dimension at least 3, evolving under the Ricci flow. The estimate depends on some 
constants arising from a Sobolev imbedding theorem. Considering the case when the scalar curvature is 
positive throughout the manifold, at any time, we will obtain, as a corollary, a bound similar to the one 
known for the fixed metric case. 

1. Introduction 

We consider a closed manifold M, whose metric evolves by the Ricci flow: 

d 

(1.1) — g(x,t) = -2Bic(x,t), xeM,te[0,T] 

and we obtain un upper bound for the heat kernel G{x, t; y, s), which is the fundamental solution of the heat 
equation 

(1.2) (A-^ju(x,t)=0, xEM,te[0,T}. 

Determining bounds for the heat operator on manifolds has been a topic of interest, as it had proven to 
have many applications. D. Aronson made use of a parabolic Harnack inequality to bound the fundamental 
solution for a general second-order parabolic operator [Aro67 . Later, in their celebrated paper (LY86] . P. 
Li and S.-T. Yau derived gradient estimates for positive solutions to the heat equation on closed manifolds 
with bounded Ricci curvature, from which they obtained Harnack inequalities. Further these inequalities 
were used to get upper and lower bounds on the heat kernel. They considered manifolds with boundaries, 
satisfying Dirichlet and Neumann boundary conditions, the heat kernel bounds extending to the boundary 
when the latter was convex. Further, J. Wang derived in |Wan97j a global version of gradient estimates when 
the boundary is nonconvex, and he obtained both upper and lower bounds for the heat kernel satisfying 
Neumann conditions. 

In geometric analysis, heat kernel estimates, together with Sobolev imbedding theorems, have been proven 
useful in the study of Ricci flows, especially in the case with surgeries. Since Sobolev imbeddings and 
inequalities relate the integrability (in some L p sense) of the derivative of a function to the integrability of 
the function itself, they become useful when looking at partial differential equations. They also have proven 
useful in characterizing the space where the function is defined (for a detailed discussion see, for example, 
[SCOT]). 

In [Gue02j . C. Guenther studied the fundamental solution of the linear parabolic operator L{u) = (A — 
— h)u, on compact n-dimensional manifolds with time dependent metric, where h is a smooth space-time 
function. She proved the uniqueness, positivity, the adjoint property and the semigroup property of this 
operator, which thus behaves like the usual heat kernel. As a particular case (h = 0), she obtained the 
existence and properties of the heat kernel under the Ricci flow. 

G. Perelman gave a proof in |Per02j of the pseudolocality theorem, which states that Euclidean looking 
regions in closed manifolds evolving by Ricci flow remain localized, under some curvature assumptions. In 
order to prove this, he obtained a differential Li-Yau-Hamilton type inequality for the fundamental solution 
of the conjugate heat equation Au + J^- — Ru = 0, where R is the scalar curvature on the manifold. Later 
S.-Y. Hsu obtained in |Hsu08l . by a variation of the method introduced by P. Li, S.-T. Yau and J. Wang, a 
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gradient estimate for the solution of the conjugate heat equation on closed manifolds under Ricci flow, and 
as a consequence, bounds for its fundamental solution. 

Q. Zhang also considered in Zha06] the conjugate heat equation introduced by Perelman, but after a time 
reversal: Au — Jj — Ru = 0, and the metric evolving under forward Ricci flow ^g(x,t) = 2Ric(x,t). He 
considered a complete manifold, with Ric(#(i)) > k and the injectivity radius i > 0. He obtained un upper 
bound on the fundamental solution of this equation using the Nash method, without any gradient estimate, 
and his result depends on the best constants in the Sobolev imbedding theorem. Let's just mention that our 
heat kernel equals the fundamental solution in Q. Zhang's paper (since G(x,t;y, s) satisfies the conjugate 
heat equation in the (y, s) variables) and our result is an improvement, since there are no conditions on the 
Ricci curvature or the injectivity radius. 

Recently, Q. Zhang and X. Cao characterized in [CZlOl the Type I singularity model of the Ricci flow by 
means of upper and lower bounds of the fundamental solution of the conjugate heat equation. 

In this paper, we obtain a bound on the heat kernel, depending on the best constants in a Sobolev 
imbedding theorem. We state the Sobolev imbedding theorems and the constants that the result will depend 
on in section ([2]), while the proof is given in section (|3|). We will conclude the paper with the special 
case, when the scalar curvature is positive at the initial time (and hence, since the Ricci flow preserves the 
positivity of the scalar curvature, throughout the manifold at all given times). 

We should add that the bound we get is not sharp, and that a long term goal is to get estimates similar 
to the Li-Yau ones. This, however, has proved to be much more difficult, due to the changing nature of the 
metric. Hope comes from the fact that similar gradient estimates on the solution of the heat equation have 
been found by the author, together with X. Cao and A. Pulemotov, using similar methods (see MCP09 ). 

The main result of this paper can be stated as follows: 

Theorem 1.1. Let M n be a closed Riemannian manifold, with n > 3 and let (M, g(x,t)^) ,t € [0,T] be a 
solution to the Ricci flow Let G(x,t;y, s) be the heat kernel, i.e. fundamental solution for the heat 

equation (| 1 - 2[) . Then there exists a positive number C n , which depends only on the dimension n of the 
manifold such that: 

C 

G(x,t;y,s) < — ; 



If \~ 2 2 »H 



(It 



for < s < t < T; here H{t) = J — | • mo * c T dr, 1/mo = inf t= o R - the infimum of the scalar 

s 

curvature, taken at time 0, and A(t) and B(t) are positive functions, which depend on the best constant in 
the Sobolev imbedding theorem. 

One can notice that there are no curvature assumptions, just like in [Zha06j and in |Hsu08j . where the 
conjugate heat equation was analysed. 

The estimate may not seem natural, but in a special case, when the scalar curvature R(x, 0) > (and 
thus R(x,t) > for any t £ [0,T]), one obtains a similar result to the one in the fixed metric case. Let's 
recall that J. Wang obtained in [Wan97j that the heat kernel on an n-dimensional compact Riemannian 
manifold M, with fixed metric, is bounded from above by C(S)(t — s) _n / 2 , where C(S) is the Neumann 
Sobolev constant of M , coming from a Sobolev imbedding theorem. Our corollary exibits a similar bound: 

Corollary 1.2. Under the same assumptions as in theorem together with the condition that the scalar 

curvature R{x, t) be (strictly) positive at t = 0, there exists a positive number C n , which depends only on the 
dimension n of the manifold and on the best constant in the Sobolev imbedding theorem in W 1 , such that: 

G(x,t;y,s)<C n - „ for < s < t < T 

(t-s)2 

The exact expression of G n will be shown in the proof. C n differs from C(S) (as it appears in |Wan97j ). 
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2. Setup 

We consider an rt-dimensional manifold without boundary M, which is compact, connected, oriented and 
smooth. 

For T > 0, let (M,g(x,t)),t 6 [0,T] to be a solution to the Ricci flow 

d 

— g(x,t) = -2Ric(x,t) 

The interval that we consider [0, T] is a subset of the interval of short-time existence, hence we won't deal 
with blow-ups. 

We say that a smooth positive function u : M x [0, T] — > R satisfies the heat equation if the following 
holds 

d \ 

A-—\u(x,t)=0 

Here, A stands for the Laplacian given by g(x,t). We will use V and | • | to denote the gradient and the 
norm with respect to g(x,t), respectively. We emphasize that A, V, and • | all depend on t G [0,T]. XY 
denotes the scalar product of the vectors X and Y with respect to the metric g(x,t). 

We will denote the heat kernel, i.e. the fundamental solution of the heat operator (A — Jj) by G(x, t; y, s). 
Let's recall that the fundamental solution of an operator L is a smooth function G(x, t; y, s) : M x [0, T] x 
M x [0, T] — > K, with s < t, which satisfies two properties: 

(i) L(G) = 0in (x,t) for (x,t) + {y,s) 

(ii) lim t _,. s G(., t; y, s) = S y for every y, where S v is the Dirac delta function. 

Guenther proved the existence and studied the properties of the fundamental solution for the operator 
L(u) = (A — — h{x,t)~) u on a compact manifold whose metric evolves under the Ricci flow (h(x,t) is a 
smooth function) |Gue02[ Theorem 2.1]. In particular, if h(x,t) = we get the existence of the heat kernel. 

During the computations, we will drop the arguments {x, t; y, s) since it will be clear with respect to which 
variables we are considering the measure over which we are integrating. 

Our proof relies on two Sobolev imbedding theorems, which are stated below. Since the manifold is 
compact, there will be no assumption on the injectivity radius or on the Ricci curvature (as in |Zha06| ). 

In |Aub76j . Aubin proved the following: 

Theorem 2.1. Let M n be a compact Riemannian manifold. If 1 < q < n, then for any e > 0, for any 
q' G [1,?] and for any r > 1, there exists a positive constant B(g,e,n) such that for any ip G W 1,q (M) (the 
Sobolev space of weakly differentiable functions) : 

IHi;< (K( n , q y + e)||v^||; + B|Mip . 

Here K(n,q) is the best constant in the Sobolev imbedding (inequality) in R™ and p = (nq) / (n — q). 

Along the Ricci flow, Zhang proved the following uniform Sobolev inequality in [Zha09 : 

Theorem 2.2. Let M n be a compact Riemannian manifold, with n > 3 and let (M, 5 , (i)) t£ r r i be a solution 
to the Ricci flow Let A and B be positive numbers such that for (M,g(0)) the following Sobolev 

inequality holds: for any v G M /1,2 (M, g(0)), 



v\*=* d[*(g(0)) <A \Vv\ z d»(9(P)) + B v d//( 5 (0)) 

M / JM JM 

Then there exist positive functions A(t), B(t) depending only on the initial metric g(0) in terms of A and 
B, and t such that, for all v G W 1,2 (M, g(t)), t > 0, the following holds 

n-2 

(J \v\& M9(t))j " < Mt) f (|V«| 2 + \Rv 2 ^j df,(g(t)) + B{t) f v 2 dfi(g(t)) 

Here R is the scalar curvature with respect to g(t). Moreover, if R(x,0) > 0, then A{t) and B(t) are 
independent oft. 
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3. Proof 



We assume without loss of generality that s = 0. By the semigroup property of the heat kernel |Gue02[ 
Theorem 2.6] and the Cauchy-Bunyakovsky-Schwarz inequality we have that: 



G(x,t;y,0) 



G ( x,t;z, ^ ) G (z,-;y,0 ) d/j, (z, ^ 



< 



M 



At 



G- { x,t;z, - j d/j,(z, ^ 



1/2 



A/ 



g- ! z, 2 ;y '°) dM ( z ' 5 



1/2 



The key of the proof consists in determining upper bounds for the following two quantities: 

a{t) = / G 2 (x,t;y,s)d[j,(x,t) (for s fixed) 
JM 



G 2 (x, t; y, s)dfi(y, s) (for t fixed) 



At 



Once we have these bounds, the conclusion follows immediately. The strategy to get these bounds consists 
in finding an ordinary differential inequality for each of the two quantities. 

First let's recall that, by definition, G satisfies the heat equation in the (x,t) coordinates 

A x G(x, t; y, s) - d t G(x, t; y,s) = 

whereas in the (y, s) it satisfies the conjugate heat equation 

A y G(x, t; y, s) + d s G(x, t; y, s) - R(y, s)G(x, t; y, s) = 

here R(y, s) is the scalar curvature, measured with respect to the metric g(s). 

We will first deduce a bound on a(t), by finding an inequality involving a'(t) and a(t). Note that we will 
treat G as being a function of (x, t), the (y, s) part is fixed. 

Since 4f{df£) — —Rdfi (due to the Ricci flow), one has: 

(3.1) a'{t) = 2 I G ■ G t dfi(x,t) - f G 2 R d(i{x, t) = 2 / G ■ (AG) dfi(x,t) - [ G 2 Rd^{x,t) = 

JM JM JM JM 

(3.2) -2 f |VG| 2 dfj,- [ G 2 Rdfi<- [ [|VG| 2 + RG 2 } dfj,(x, t) 

JM JM JM 

The difficult part will be to estimate J M \VG\ 2 d/j,. The way we proceed is to use the Sobolev imbedding 
theorem, which gives a relation between J M |VG| 2 d/x and J M G 2 d[i, and the Holder inequality to bound the 
term involving Q 2n /( n - 2 ) : 



(3.3) 



G 2 dfx(x,t) < 



A I 



G™- 2 dfj,(x,t) 



At 



At 



G dfi(x, t) 



By theorem (|2.1[) , for r = 2, q = q 1 = 2 and p = 2n/(n — 2) one gets that at time t = 0, the following 
inequality holds for any v £ W 1,2 (M t g(0)) and for some B > 0: 



\v\^ dn(g(0))) " <(K(n,2) 2 + e) \ Vv\ 2 dfi(g(0)) + B / v 2 dfi(g(0)) 

lAI J JM JM 

where K(n, 2) - the best constant in the Sobolev imbedding in R. 

From this, by theorem 12.21 it follows that at any time t > within the interval of short existance of the 
solution to the Ricci flow [0,T) one has that for all v 6 W 1,2 {M, g(t)): 



\v\&* dfi(g(t)) ) " <A(t) f ( \Vv\ 2 + -Rv 2 ) dfi(g(t))+B(t) [ v 2 dp(g(t)) 

IM / J At \ 4 / J At 

where A(t) is a positive function depending on g(0) and (K(n, 2) 2 + e), while B(t) is also a positive 
function, depending on B. 
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Since G(x,t; ., .) e H /1,2 (M, g(t)) (it is even smooth), then the above holds, so one can relate the RHS of 
(|3.3[) to the Sobolev inequality: 

(3.4) 



G 2 dfi(x,t) < 



M 



G™" 2 dfi(x,t) 



M 



G dfx(x, t) 



M 



(3.5) 



< 



A(t) I |VG| 2 + -RG 2 dfj,{x,t) + B(t) / G dp(x,t) 
m V 4 / J M 



G d[x(x, t) 



M 



We need to estimate the term J(t) := J M G(x,t;y, s) dfJi(x,t). By the definition of the fundamental 
solution, we have that: J M G(x,t;y, s)dfj,(y, s) = 1, but that's not true if one integrates in (x,t). We will 
obtain a differential inequality for J(t) and the estimate will follow therefrom. 

J'(t)= I G t (x,t;y,s) dfi(x,t) + j G(x,t;y,s)j-^ dfj,(x,t) = / A x G(x,t;y,s) dfi(x,t) 



M 



M 



G(x,t;y, s)R(x,t) dfi(x,t) = — / G(x,t;y, s)R(x,t) dfi(x,t) 

M JM 

the first term being 0, as we are on a compact manifold, without boundary. 

The scalar curvature satisfies the following differential inequality (see CLN06J): 

ot n 

Since the solutions of the ODE = ^p 2 are p(t) — np ^-i _ 2t ; by the maximum principle we get a bound 
on the scalar curvature, for s < r < t: 



R{z,t) > 



1 



n(inf t=0 R)- 1 - 2t (inf i=0 i?)' 1 - }r 
(here and later, if inf t=0 R > 0, then the above is regarded as zero). 

Using this lower bound for R (for r £ (s,t]), we get: 



TO - C n T 



1 



-J(r) 



J'(r) < 

m ~ c n T 

After integrating the above from s to t, while noting that by J(s) one understands: 

J(s) = lim / G(x,t;y,s) dfi(x,t) — / limG(x,t;y, s) dfi(x, t) = / 5 y (x) dp(x, s) = 1 



A I 



M ' 



M 



one obtains: 



J(t)<(^) 7 H^ 

\mo- c n s J 

Hence J M G(x,t;y, s) dfi(x,t) < (xt.s))^ and (|3.4[) becomes: 



G 2 dn{x,t) < 



M 



A(t) j \^\VG\ 2 + -RG 2 ) dp{x,t) + B{t) / G d[i{x,t) 



M 



{xt.sY 



From this it follows immediately that: 

1 



|VG| 2 dfi(x,t) > 



M 



xLA(t) 



G 2 dfj,(x,t) 



B(t) 



G 2 dn(x,t) - - / RG 2 d/nr.1) 



M j A(t) J M 

Combining this with the inequality from (|3.ip . one obtains the following differential inequality for a(t): 

1 



n ' j/j ' " X 2 t A(t) a{t) ^ + W) ait) ~ i ' /lY72 ' / /' , - , '- /) 
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Note that the above is true for any r G (s,t\. For the following computation, we will consider t fixed as 



well. Recall that for r <E [s,t], R(-,t) > 



h(r) :- 



Denoting with: 



B(t) 3 



1 



A(t) 4 m - c„t 



we get: 



«'(r) < - 1 a{T)^ +h{r)a{T) 



Let -ff(r) be an antiderivative of ft-(r). By the integrating factor method, one finds: 



{e- H ^a{ T ))' < — 



1 



X 2 {t)A(t) 

Since the above is true for any r G (s,t\, by integrating from s to t and taking into account that 
lim a(r) = / lim G 2 (x, r; y, s) dfj,(x, t) = / 5 2 {x) dfi(x, s) = 



one obtains the first necessary bound: 



a(t) < 



I ^(r)A(r) dT 



where C n = (|) 



1U 



The next step is to estimate /3(s) = J M G 2 (x, t; y, s) d/i(y, s), for which the computation is different, due 
to the assymetry of the equation. As stated above, the second entries of G satisfy the conjugated equation: 

A y G(x, t; y, s) + d s G(x, t; y, s) - RG(x, t; y, s) = 

Proceeding just as in the a(t) case, we get the following: 

P'(s) = 2[ GG S d/j,{y,s) - [ RG 2 d/j,{y , s) = 2 f G(-AG + RG) dfj,{y,s) - [ RG 2 dfi(y,s) 
Jm Jm Jm Jm 

= -2 [ G(AG) dn(y,s)+ f RG 2 dfi(y, s) = 2 f |VG| 2 dfi(y, s) + [ RG 2 dfi(y, s) 
Jm Jm Jm Jm 

> f |VG| 2 dn(y,s)+ f RG 2 d^{y,s) 

JM JM 

Hence 

P'(s)> I (\VG\ 2 + RG 2 )d^y,s) 

JM 

But this time, by the property of the heat kernel: 

J(s) := / G(x, t; y, s) dfi(y, s) = 1 

JM 

so by applying Holder (as for a{t)) and relating it to the Sobolev inequality, we get: 



JM 



G 2 dn(y, s) < 



A(s) 



M 



|VG| 2 + -RG 2 ) dfi(y, s) + B(s) I G 2 dn(y, s) 



M 



M 



G d/i(y, s) 



A(s) J (j VG| 2 + - A RG 2 ^j dn(y, s) + B(s) J G 2 d»(y, 



s) 



Following the same steps as for a(t), one finds 



ns)>^(s)^-h{sm8) 



(h(s) denotes, as before, - § • ^p^) 



BOUNDS ON THE HEAT KERNEL UNDER THE RICCI FLOW 



7 



The above is true for any t G [s,t). We will apply again the integrating factor method, with H(t) being 
the same antiderivative of h(r) as above. For re [s,t), the following holds: 

(e ffM /3(r))'> 7 i-(e^/3(r))^ e -*^) 
A(t) 

Integrating between s and t, and taking into account that 

lim/3(r)=/ lim G 2 (x, t- y, t) d^(y,T)= [ 6 2 y {x) dfi(y, t) = 



we get the second desired bound: 



P(s)< 



From the estimates of a and (3 we get the following: 



n/2 



"'I 



J^G 2 (z,l;y,<?) dfi(z,t/2)< 



't/2 

r I mg-c„T \ 
J \ m J 



dr 



C,e-« 



n/2 



i J A(t) 
V/2 



dr 



t/2 



B(t) _ 3 



A(t) 4 rriQ — CnT 



Here, we may choose H(t/2) = J 

v ' A(r) 4 mQ~c n T 

The conclusion follows from multiplying the relations above 



dr, since the relation is true for any antiderivative 



4. Special case: positive scalar curvature 

In the special case when R(x,t) > 0, one gets that J'(r) < 0, which means that J(r) is decreasing, so 
J(t) < J(s) = 1, thus leading to the differential inequality for a(t) to be: 

a ^^~W) a{t) ^ + W) a[t) 

And from this the bound for a{t) becomes: 

C n e H W 



a{t) < 



where H(t) is the antiderivative of such that H(s) ^ and H(t) ^ 0. 
Similarly, one obtaines for j3(s): 



and from this: 
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000 < 



C n e 



I 



-H(s) 



dr 



i/2 



where H(t) is the same antiderivative of -^fy as above. 

By (|2.2|) . in the case of R(x, 0) > 0, the two functions A(t) and B{t) are constants, let's call them A and 
-B. Recall that A is in fact K(n, 2) + e, where K(n, 2) is the best constant in the Sobolev imbedding and 
e > 0. 



One has that H (t) = ^t. Using this, we get: 



Cj7 



G(x,t;y,s) < 



But, by Taylor expansion, the last expression is bounded by 




2B t^j 
e rtA 2 



Cr, 



■K (i 



e tia 2 



where C n = C n ■ (2A)* = (f ) 2 • (2(tf(n,2) + e))*. 
Combining the two, we ge the desired corollary. 

References 

[Aro67] D. G. Aronson. Bound for the fundamental solution of a parabolic equation. Bull. Am. Math. Soc, 73:890-896, 1967. 
[Aub76] T. Aubin. Problemes isoperimetriques et espaces de sobolev. J. Differential Geometry, 11:573-598, 1976. 
[CLN06] B. Chow, P. Lu, and L. Ni. Hamilton's Ricci flow. American Mathematical Society, Providence, RI; Science Press, 
New York, 2006. 

[CZ10] X. Cao and Q. Zhang. The conjugate heat equation and ancient solutions of the ricci flow. arXiv:math/1006.0540vl, 
2010. 

[Gue02] CM. Guenther. The fundamental solution on manifolds with time-dependent metrics. J. Geom. Anal., 12:425-436, 
2002. 

[HsuOS] S.-Y. Hsu. Some results for the perelman lyh-type inequality. 2008. 

[LY86] P. Li and S.-T. Yau. On the parabolic kernel of the schroedinger operator. Acta Math., 156:153—201, 1986. 
[MCP09] M.Bailesteanu, X. Cao, and A. Pulemotov. Gradient estimates for the heat equation under the ricci flow. J. Funct. 
Anal, 2009. 

[Per02] G. Perelman. The entropy formula for the ricci flow and its geometric applications. arXiv:math/0211159vl, 2002. 
[SC01] L. Saloff-Coste. Aspects of Sobolev-Type Inequalities. London Mathematical Society Lecture Note Series; Cambridge 

University Press,, 2001. 
[Wan97] J. Wang. Global heat kernel estimates. Pacific J. Math., 178:377-398, 1997. 

[Zha06] Q. Zhang. Some gradient estimates for the heat equation on domains and for an equation by perelman. Int. Math. 
Res. Not., 2006:1-39, 2006. 

[Zha09] Q. Zhang. Sobolev inequalities, heat kernel estimates and applications to Ricci flow, book proposal, 2009. 



